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ABSTRACT 

We examine closely the solar Center-to-Limb variation of continua and lines and compare 
observations with predictions from both a 3-D hydrodynamic simulation of the solar surface 
(provided by M. Asplund and collaborators) and 1-D model atmospheres. Intensities from the 
3-D time series are derived by means of the new synthesis code Asset, which overcomes limita- 
tions of previously available codes by including a consistent treatment of scattering and allowing 
for arbitrarily complex line and continuum opacities. In the continuum, we find very similar 
discrepancies between synthesis and observation for both types of model atmospheres. This is 
in contrast to previous studies that used a "horizontally" and time averaged representation of 
the 3-D model and found a significantly larger disagreement with observations. The presence 
of temperature and velocity fields in the 3-D simulation provides a significant advantage when 
it comes to reproduce solar spectral line shapes. Nonetheless, a comparison of observed and 
synthetic equivalent widths reveals that the 3-D model also predicts more uniform abundances 
as a function of position angle on the disk. We conclude that the 3-D simulation provides not 
only a more realistic description of the gas dynamics, but, despite its simplified treatment of 
the radiation transport, it also predicts reasonably well the observed Center-to-Limb variation, 
which is indicative of a thermal structure free from significant systematic errors. 
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A few years back, it was realized that one of 
the most 'trusted' absorption lines to gauge the 
oxygen abundance in the solar photosphere, the 
forbidden [01] line at A6300, was blended with 
a Ni I transition. These two transitions overlap 
so closely that only a minor distortion is appar- 
ent in the observed feature. Disentangling the 
two contributions with the help of a 3-D hydro- 
dynamical simulation of surface convection led us 
to propose a reductio n of the solar photospheri c 
abundance by ~ 30% (| Allende Prieto et al.ll200lh . 
Using the same solar model, subsequent analysis 
of other atomic oxygen and OH lines confirmed 



the lower abundance, re sulting in an average value 
log <@ (0)= 8.66 ± 0.05 l| Asplund et al.ll2004 b 

This reduction in the solar O/H ratio, together 
with a parallel downward revision for carbon 



(j Allende Prieto et al.ll2002H Asplund et al.ll2005bl ). 
ruins the nearly perfect agreement between models 
of the solar int er ior an d seismological observations 
(Bah call et al. 2005i iDelahave fc Pinsonneaultl 
20061 : iLin et aTll2007l h A br ief overview of t he pro - 
posed solutions is given by I Allende Prieto ( 2007 ). 
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Interior and surface models appear to describe two 
different stars. 

Supporters of the new hydrodynamical models, 
and the revised surface abundances, focus on their 
strengths: they include more realistic physics, and 
are able to reproduce extremely well detailed ob- 
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servations (oscillations, spectral line asymmetries 
and net wavelength blueshifts, granulation con- 
trast and topology). Detractors emphasize the 
fact that the new models necessarily employ a sim- 
plified description of the radiation field and they 
have not been tested to the same extent as clas- 
sical 1-D models. The calculation of spectra for 
3-D time-dependent models is a demanding task, 
which is likely the main reason why some funda- 
mental tests have not yet been performed for the 
new models. 

On the basis of 1-D ra diative transfer calcula- 
tions, Avres et al. ( 20061 ) suggest that the ther- 
mal profile of the solar surface convection simula- 
tion of Asp l und et al. (2000) may be incorrect. 
Avres et al. ( 20061) make use of a 1-D average, 
both 'horizontal' and over time, of the 3-D simula- 
tion to analyze the center-to-limb variation in the 
continuum, finding that the averaged model per- 
forms much more poorly than the semi-empirical 
FAL C model of Fontenla et al. (1993). When the 
FAL C model is adopted, an analysis of CO lines 
leads to a much higher oxy gen abundance, and 



therefore lAyres et al.l (|2006l ) question the down 



ward revision proposed earlier. 



Asplund et al.l (|2005bl ) argue that when clas- 



sical 1-D model atmospheres are employed, the 
inferred oxygen abundance from atomic features 
differs by only 0.05 dex between an analysis in 
1-D and 3-D. The difference is even smaller for 
atomic carbon lines. When the hydrodynamical 
model is considered, there is good agreement be- 
tween the oxygen abundance inferr ed from atomic 



lines and from O H transitions ([Asplund et al 
2004t IScott et al.l l2006h . A high value of the 



oxygen abundance is derived only when consid- 
ering molecular tracers in one dimensional at- 
mospheres, perhaps not a surprising result given 
the high sensitivity to temperature of the molec- 
ular dissociation. A low oxygen abundance 
(loge(O) = 8.63) value is also deduced from 
atomic lines and atmospheric models based on the 
i nversion of spatially resolved polarimetric data 
(|Socas- Navarro fc Nortor]l2007t ). 

Despite the balance seems favorable to the 3- 
D models and the low values of the oxygen and 
carbon abundances, a failure of the 3-D model to 
match the observed li mb darkening as su ggested 
by the experiments of Avres et al. ( 20061 ). would 
be reason for serious concern. In the present pa- 



per, we perform spectral synthesis on the solar 
surface convection simulation of Asplund et al. 
(2000) with the goal of testing its ability to re- 
produce the observed center-to-limb variations of 
both the continuum intensities, and the equivalent 
widths of spectral lines. We compare its perfor- 
mance with commonly-used theoretical 1-D model 
atmospheres. Our calculations are rigorous: they 
take into account the four-dimensionality of the 
hydrodynamical simulation: its 3-D geometry and 
time dependency. After a concise description of 
our calculations in Section 2, §3 outlines the com- 
parison with solar observations and §4 summarizes 
our conclusions. 

2. Models and Spectrum Synthesis 

We investigate the Center-to-Limb Variation 
(CLV) of the solar spectrum for the continuum and 
lines. Snapshots taken from 3-D h ydrodynamical 



simula tions of the solar surface by lAsplund et al 



(2000) serve as model atmospheres. The syn- 
thetic continuum intensities and line profiles are 
calculated by means of the new spectrum syn- 
thesis code AsSeT (Advanced Spectrum Synthe- 
sis 3-D Tool), which is designed to solve accu- 
rately the equation of radiation transfer in 3-D. 
Th e new synthesis code will be described in detail 



by iKoesterke et al.~ ( in prep. ) and only the key 



features are highlighted in subsequent sections. 

2.1. Hydrodynamic Models 

The simulation of solar granulation was carried 
out with a 3-D, time-depend ent, compressible, 
radiative-hydrodynamics cod e dNordlund fc Steinl 
ll990HStein fc Nordlundlll989l : lAsplund et al.ll999l ). 
The simulation describes a volume of 6.0 x 6.0 x 3.8 Mm 
(about 1 Mm being above r cont « 1) with 200 x 200 x 82 
equidistantly spaced grid points over two hours of 
solar time. About 10 granules are included in the 
computed domain at any given time. 

99 snapshots were taken in 30 s intervals from 
a shorter sequence of 50 min. The grid points and 
the physical dimensions are changed to accommo- 
date the spectrum synthesis: The horizontal sam- 
pling is reduced by omitting 3 out of 4 grid points 
in both directions; the vertical extension is de- 
creased by omitting layers below t™JJJ w 300 while 
keeping the number of grid points in z-direction 
constant, i.e. by increasing the vertical sampling 
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and introducing a non-equidistant vertical grid. 
After these changes, a single snapshot covers ap- 
proximately a volume of 6.0x6.0xl.7Mm with 
50x50x82 grid points (lAsplund et alj|2000h . 



2.2. Spectrum Synthesis 

Compared to the spectrum synthesis in one di- 
mension, the calculation of emergent fluxes and in- 
tensities from 3-D snapshots is a tremendous task, 
even w hen LTE is app li ed. Pr ev ious investigations 



e.g., lAsplund et al 



(l2000h : iLudwig fc Steffenl 
(|20Q7l ) ) were limited to the calculation of a single 
line profile or a blend of very few individual lines 
on top of constant background opacities, and with- 
out scattering. In order to overcome these limita- 
tions, we devise a new scheme that is capable of 
dealing with arbitrary line blends, frequency de- 
pendent continuum opacities, and scattering. The 
spectrum synthesis is divided into five separate 
tasks that are outlined below. A more detailed 
description which co ntains all essential numerica l 
tests will be given bv lKoesterke et al~ ( in prep. ). 



2.2.1. Opacity Interpolation 

For the 3D calculations we face a situation in 
which we have to provide detailed opacities for 
ss 2-10 7 grid points for every single frequency under 
consideration. Under the assumption of LTE, the 
size of the problem can be reduced substantially 
by using an interpolation scheme to derive opac- 
ities from a dataset that has orders-of-magnitude 
fewer datapoints. We introduce an opacity grid 
that covers all grid points of the snapshots in the 
temperature-density plane. The grid points are 
regularly spaced in log T and log p with typical in- 
tervals of 0.018 dex and 0.25 dex, respectively. 

We use piecewise cubic Bezier poly nomials tha t 
do not introduce artificial extrema ( Auerl [20031 . 
To enable 3 rd -order interpolations close to the 
edges, additional points are added to the opacity 
grid. The estimated interpolation error is well be- 
low 0.1% for the setup used throughout the present 
paper. 

2.2.2. Opacity Calculation 

We use a modified version of SYNSPEC 
dHubenv fc Land Il995h to prepare frequency- 
dependent opacities for the relatively small num- 
bers of grid points in the opacity grid. The mod- 



ifications allow for the calculation of opacities on 
equidistant log(A) scales, to output the opacities 
to binary files, and to skip the calculation of in- 
tensities. 

Two datasets are produced. Continuum opac- 
ities are calculated at intervals of about 1 A at 
3000 A. Full opacities (continuum and lines) are 
provided at a much finer spacing of 0.3t> m i n , with 
fmin being the thermal velocity of an iron atom at 
the minimum temperature of all grid points in all 
snapshots under consideration. A typical step in 
wavelength is 2.7 • 10" 3 A at 3000 A, which corre- 
sponds to 0.27 kms -1 . 

We adopt the sola r photospheric a bundan ces re- 
cently proposed by lAsplund et al.1 ( 2005al ). with 
carbon and oxygen abundances of loge = 8.39 
and 8.66, respectively, which are about 30% less 
than in earlier compilations ( Grevesse fc Sauval 
19981 ). We account for bound-bound and bound- 
free opacities from hydrogen and from the first two 
ionization stages of He, C, N, O, Na, Mg, Al, Si, 
Ca and Fe. Bound-free cross sections for all metals 
but iron are tak en from TOPBASE and smo othed 
as described bv lAllende Prieto et al~ ( 2003h . Iron 
bound-free opacities are derived from the pho- 
toionization cross -sections computed by t he Iro n 
Project (see, e.g jNaharl (|l995l ); iBautistal (|l997l )l. 
after smoothing. 

Bound-bound log (gf) values are taken from 
Kurucz, augmented by damping constants from 
Barklem et al. I (l2000h where available. We also ac- 
count for bound- free opacities from H~, H^, CH 
and OH, and for a few million molecular lines from 
the nine most prominent molecules in the wave- 
length range from 2200 A to 7200 A. Thomson and 
Rayleigh (atomic hydrogen) scattering are consid- 
ered as well, as described below in £)2.2.31 The 
equation of state is solved considering the first 
99 elements and 338 molecular species. Chemi- 
cal equilibrium is assumed for the calculation of 
the molecular abundances, and the atomic abun- 
dances are updated accordingly (private comm. 
from I. Hubeny). 

2.2.3. Scattering 

We employ a background approximation, cal- 
culating the radiation field J v for the sparse con- 
tinuum frequency points for which we have cal- 
culated the continuum opacity without any con- 
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tribution from spectral lines. The calculation 
starts at the bluemost frequency and the velocity 
field is neglected at this point: no frequency cou- 
pling is present. The opacities for individual grid 
points are derived by interpolation from the opac- 
ity grid, and the emissivities are calculated assum- 
ing LTE. As mentioned above, we include electron 
(Thomson) scattering and Raylcigh scattering by 
atomic hydrogen. An Accelerated Lambda Itera- 
tion (ALI) scheme is used to obtain a consistent 
solution of the mean radiation field J v and the 
source function S v at all grid points. In turn, J v is 
calculated from S v and vice versa, accelerating the 
iteration by amplifying AJ U = J^ cw - J° ld by the 
factor 1 / (1 — A*) with A* being the app roximate 
lambda operator ( Olson fc Kunasd 1987 ). Gener- 
ally the mean radiation from the last frequency 
point, i.e. the frequency to the blue, serves as 
an initial guess of J v at the actual frequency. At 
the first frequency point, the iteration starts with 

The formal solution, i.e. the solution of the 
equation J v = AS^, is obtai ned by means of 
a sho rt characteristics scheme ( Olson k, Kunaszl 
19871 ). For all grid points the angle-dependent 
intensity Ijf is derived by integrating the source 
function along the ray between the grid point it- 
self and the closest intersection of the ray with a 
horizontal or vertical plane in the mesh. The op- 
erator A* needed for the acceleration is calculated 
within the formal solution. 

For the present calculations, J v is integrated 
from I£ at 48 angles (6 in /i, 8 in 0). The integra- 
tion in \x is performed by a three-point Gaussian 
quadrature for each half-space, i.e. for rays point- 
ing to the outer and the inner boundary, respec- 
tively. The integration in <f> is trapezoidal. The 
opacities and source functions are assumed to vary 
linearly (l st -order scheme) along the ray. 

In order to integrate the intensity between the 
grid point and the point of intersection where the 
ray leaves the grid cell, the opacity, source func- 
tion and the specific intensity (/€„, S v , I£) have 
to be provided at both ends of the ray. Since 
the point where the ray leaves the cell is generally 
not a grid point itself, an interpolation scheme has 
to be employed to derive the required quantities. 
We perform interpolations in two dimensions on 
the surfaces of the cuboids applying again Bezier 
polynomials with control values that avoid the in- 



troduction of artificial extrema. The interpolation 
may introduce a noticeable source of numerical in- 
accuracies. Detailed tests, using an artificial 3-D 
structure constructed by horizontally replicating 
a 1-D model, revealed that a 3 rd -order interpola- 
tion scheme provides sufficient accuracy where lin- 
ear interpolations fail in reproducing the radiation 
field: the mean relative errors are 0.5% and 0.05% 
for linear and cubic interpolation, respectively. 

It is possible, in terms of computing time, to 
calculate J v from the full opacity dataset for all 
frequencies (our 'fine' sampling). However, since 
the total effect of scattering for the solar case in 
the optical is quite small, the differences between 
the two methods are negligible. Therefore, we ap- 
ply the faster method throughout this paper. Note 
that in both approximations (using background or 
full opacities), the calculation of the mean radia- 
tion field does not account for any frequency cou- 
pling. 

2.2.4- Calculation of Intensities and Fluxes 

The emergent flux is calculated from the opac- 
ities of the full dataset provided at the fine fre- 
quency grid. Again, the opacities for individual 
grid points are derived by interpolation from the 
opacity grid and the emissivities are calculated 
from LTE. The mean background radiation field 
J„ is interpolated from the coarser continuum fre- 
quency grid to the actual frequency, and it con- 
tributes to the source function at all grid points via 
Thomson and Rayleigh (atomic hydrogen) scatter- 
ing opacities. 

The integration along a ray is performed in 
the observer's frame by following long character- 
istics from the top layer down to optical depths 
of TR, ay > 20. Frequency shifts due to the veloc- 
ity field are applied to the opacities and source 
functions. Each ray starts at a grid point of the 
top layer and is built by the points of intersec- 
tion of the ray and the mesh. At these points of 
intersection an interpolation in three dimensions 
is generally performed, i.e. a 2-D geometric in- 
terpolation in the X-Y, X-Z or Y-Z-plane, respec- 
tively, is enhanced by an interpolation in frequency 
necessitated by the presence of the velocity field. 
Additional points are inserted into the ray to en- 
sure full frequency coverage of the opacities. This 
is done when the difference of the velocity field 
projected onto the ray between the entry and exit 
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point of a grid cell exceeds the frequency spacing 
of the opacity. Without these additional points 
and in the presence of large velocity gradients, 
line opacities could be underestimated along the 
ray -a line could be shifted to one side at the en- 
try point and to the other side at the exit point-, 
leaving only neighboring continuum opacities vis- 
ible to both points while the line is hidden within 
the cell. 

Similar to the calculation of the mean radiation 
field J described in Sect. 12.2.31 all interpolations 
in both space and frequency are based on piece- 
wise cubic Bezier polynomials. It is not completely 
trivial to mention that for the accurate calculation 
of the emergent intensities, the application of a 
high-order interpolation scheme is much more im- 
portant than it is for the calculation of the mean 
background radiation field (Sect. I2.2.3[) . Here we 
are calculating precisely the quantity we are inter- 
ested in, i.e. specific intensities. But, in addition 
to that, we deal with interpolations in three di- 
mensions (2-D in space, 1-D in frequency) instead 
of a 2-D interpolation in space. Hence, any quan- 
tity is derived from 21 1-D interpolations rather 
than just 5. 

In the standard setup of the 3-D calculations, 
20 rays are used for the integration of the flux F v 
from the intensities Ijf. Similar to the integration 
of J described in Sect. I2.2.3) the integration in 

is a three-point Gaussian quadrature, while the 
integration in (j> is trapezoidal. Eight angles in <f> 
are assigned to the first two of the fi angles while 
the last and most inclined angle with the by far 
smallest (flux) integration weight has 4 contribut- 
ing cj> angles. Note that for the investigation of 
the Center-to-Limb variation, the number of an- 
gles and their distribution in fj, and <f> differs con- 
siderably from this standard setup, as explained 
below (Sect. E]). 



2.3. Spectra in 1-D 

To facilitate consistent comparisons of spectra 
from 3-D and 1-D models, the new spectrum syn- 
thesis code AsseT accepts also 1-D structures as 
input. Consistency is achieved by the use of the 
same opacity data (cf. Sect. 12.2.2)) and its inter- 
polation (if desired in 1-D, cf. Sect. 12.2. ip and 
by the application of the same radiation transfer 
solvers, i.e. l st -order short and long-characteristic 
schemes (cf. Sect. [2~273l and l2.2.4[ respectively). 



All angle integrations are performed by means of 
a three-point Gaussian formula. This leaves the 
interpolations inherent to the radiation transfer 
scheme in 3-D as the only major inconsistency be- 
tween the spectra in 1-D and 3-D. Numerical tests 
have revealed that these remaining inconsistencies 
are quite small, as we will report in an upcoming 
paper. 

2.4. Solar Model 

Our choice is not to use a semi-empirical model 
of the solar atmosphere as a 1-D comparison with 
the 3-D hydrodynamical simulation, but a theo- 
retical model atmosphere. Semi-empirical models 
take advantage of observations to constrain the 
atmospheric structure, a fact that would consti- 
tute an unfair advantage over the 3-D simulation. 
Some semi-empirical models, in particular, use ob- 
served limb darkening curves, and of course it is 
meaningless to test their ability to reproduce the 
same or different observations of the center-to- 
limb variation in the continuum. Consequently 
we are using models from Kurucz, the MARCS 
group, and a horizontal- and time-averaged repre- 
sentation of the 3-D hydrodynamical simulations. 

We have derived a 1-D solar refe rence model 
from the Kurucz grid ( Kuruczl 1993f ). The refer- 
ence model is derived from 3 rd -order interpolations 
in t, T e ff, logg, and Z. Details of the interpo- 
lation scheme will be presented elsewhere. We 
have adopted the usual values of T c g — 5777 K 
and logg = 4.437 (cgs) but a reduced metallic- 
ity of log(Z/Z Q ) = —0.2 in an attempt to ac- 
count globally for the difference between the so- 
lar abundances (mainly iron) used in the calcula- 
tion of the model and more rec e nt va lues, as de- 
scribed bv lAllende Prieto et al. ( 20061 ). To avoid 
a biased result by using a single 1-D comparison 
model, we have also experimented with a solar 
MARCS model kindly provided by M. Asplund, 
and a solar model interpolated from the more re- 
cent ODFNEW grid from Kurucz (available on his 
websited). No metallicity correction was applied 
to these newer solar models. 

In earlier investi gations (|Avres et al"1l200(j|) . a 
1-D representation (| Asplund et al.ll2005bl ) of the 
3-D time series, i.e. a 'horizontal' averag^l over 



^http: / /kurucz. harvard.edu/ 

4 Horizontal average, in this context, refers to the mean 
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time, has been used to study the thermal profile 
of the 3-D model. While this approximation allows 
easy handling by means of a 1-D radiation transfer 
code, the validity of this approach has never been 
established. In order to investigate the limitations 
of this shortcut, we compare its Center-to-Limb 
variation in the continuum with the exact result 
from the 3-D radiation transfer on the full scries 
of snapshots. 

3. Center-to-Limb Variation 
3.1. Continuum 



Neckel fc Labi (11994 : iNeckell (|2003l . 120051) in- 
vestigated the Center-to-Limb variation of the 
Sun based on observations taken at the Na- 
tional Solar Observatory/Kitt Peak in 1986 and 
1987. They describe the observed intensities 
across the solar disk as a function of the he- 
liocentric distance by 5 th -order polynomials for 
30 frequencies between 303 nm an d 109 9 nm. 
Similar observa ti ons b y Petro et al. (1984) and 
Elste fc Gilliam! ( 2007) (with a smal l er spe ctral 
coverage) indicate that Neckel k, Labd ( 1994 ) may 
have over corrected for scattered light, but con- 
firm a level of accuracy of wO.4%. We have 
calculated fluxes and intensities for small spec- 
tral regions (±5kms _1 ) around eight frequen- 
cies (corresponding to standard air wavelengths 
of 3033.27A, 3499.47A, 4163.19A, 4774.27A, 
5798.80A, 7487.10A, 8117.60A, 8696.00A) and 
compare mon ochromatic synthetic in tensities with 
the data from Neckel &: Labsl ( 1994). Because the 
spectral regions are essentially free from absorp- 
tion lines, the width of the bandpass of the ob- 
servations varying between 1.5 kms^ 1 in the blue 
(3030A) and 1.9 kms" 1 in the red (10990A) is 
irrelevant. 

The fluxes were integrated from 20/3 angles (cf. 
Sects. I2T2T41 and |23|) for the 3-D/l-D calculations, 
respectively. For the study of the CLV, intensities 
(as a function of fx) were calculated for 11 posi- 
tions on the Sun (/j = cos0 = 1.0, 0.9, ...,0.1, 0.05) 
averaging over 4 directions in <j) and all horizontal 
(X-Y) positions. All 99 snapshots were utilized for 
the 3-D calculations. 

The eight frequencies cover a broad spectral 



range. Although some neighboring features are 
poorly matched by our synthetic spectr a, the so- 
lar flux spectrum of iKurucz et~ai] (|l984h is repro- 
duced well at the frequencies selected by Neckel 
& Labs, and therefore modifications of our linelist 
were deemed unnecessary (see Fig. [lj . The nor- 
malization of the synthetic spectra was achieved 
by means of "pure-continuum" fluxes that were 
derived from calculations lacking all atomic and 
molecular line opacities - Fig. [T] shows that Neckel 
& Labs did a superb job selecting continuum win- 
dows. 

Comparisons of observed and synthetic CLV's 
are conducted with datasets that are normalized 
with respect to the intensity at the disk center, i.e. 
all intensities are divided by the central intensity. 
We show the residual CLV's 



R-CLV 



value over a surface with constant vertical optical depth, 
rather than over a constant geometrical depth. 



if/^-ijyv/q?!, (i) 

in Fig. O The R-CLV's within each group are 
quite homogeneous. In addition to the data de- 
rived from our 1-D Kurucz model (cf. Sect. 12. 4| . 
we show also data from two other 1-D models, 
i.e. a MARCS model (1 st panel) from Asplund 
(priv. comm.) and an alternative (odfnew) Ku- 
rucz model (2 nd panel) from a different model 
grid (http://kurucz.harvard.edu/grids.html). The 
Center-to-Limb variation from both alternatives 
show much larger residuals and are not used for 
the comparison with the 3-D data. However, the 
scatter within the 1-D data demonstrates vividly 
the divergence that still persist among different 
1-D models. 

Our reference 1-D model (3 rd panel) describes 
the observed CLV's well down to /i « 0.5. Closer 
to the rim the R-CLV's rise to w 0.1 at /i = 0.2 
followed by a sharp decline at the rim. In 3-D (4 th 
panel) we find on average a linear trend of the R- 
CLV's with n, showing a maximum residual of > 
0.2 close to the rim. 

The investigation of the Center-to-Limb varia- 
tion of the continuum is an effective tool to probe 
the continuum forming region at and above r « 1 . 
Deviations from the observed CLV's indicate that 
the temperature gradient around r « 2/3 is incor- 
rectly reproduced by the model atmosphere. This 
can, of course, mean that the gradient in the model 
is inaccurate, but it can also signal that the opac- 
ity used for the construction of the model atmo- 
sphere differs significantly from the opacity used 
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Fig. 1. — Comparison of the normalized solar spectrum (thin, solid) and the synthetic spectra from the 1-D 
Kurucz model (dashed) and the 3-D Hydro-simulation (thick, solid) for the 8 wavelengths under consider- 
ation. The 3-D calculations were performed only for the small windows of ±5 km s _1 used here to study 
the center-to-limb variation in the continuum. For the normalization, the synthetic spectra were divided 
by the corresponding "pure-continuum" spectra. The 1-D spectrum has been convolved with a Gaussian of 
FWHM= 4.3 km s _1 , to account for macro-turbulence (FWHM= 4.2 km s _1 ) and the finite resolution of 
the solar atlas (FWHM= 0.8 km s -1 ). Note that the log-gf values of the individual lines have not been 
adjusted to match the spectrum. 
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Fig. 2. — Residual CLV's (R-CLV's) in the continuum, i.e. the difference of the observed and the synthetic 
normalized CLV's for the 3-D (solid) and the 1-D (dashed) model. The upper four panels show the 1-D 
and 3-D data combined, respectively. Average values are indicated by circles. The lower panels compare 
separately the data for the 6 wavelengths under consideration but do not repeat the data from the alternative 
(ODFNEW) Kurucz and the Marcs model. Positive R-CLV's indicate that the temperature drops off too 
fast in the model atmospheres. 
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for the spectrum synthesis. In that case, the tem- 
perature gradient is tested at the wrong depth due 
to the shift of the r-scales. 

Our spectrum calculations suffer from an incon- 
sistency introduced by the fact that the abundance 
pattern and the opacity cross-sections might dif- 
fer from what was used when the model was con- 
structed. In our reference 1-D model, we com- 
pensate for the new solar iron abundance (ep c : 
7.63 -► 7.45) and interpolate to log (Z/Z e ) = 
-0.2 in the Kurucz model grid (cf. Sect. 12. 4j) . 
The 3-D model has been c onstructed based on the 
Grevesse fe Sauyal (Il998 ) solar abundances (cf. 
Asplund et all (|2000h ) with e Fc = 7.50 and, to 
first order, no compensation is necessary. (And 
the same is true for the other two 1-D models con- 
sidered in Fig. [21) The changes in carbon and 
oxygen abundances do not affect the continuum 
opacities, which are dominated in the optical by 
H and H _ . Consequently, only metals that con- 
tribute to the electron density and therefore to the 
H~ population (i.e. Fe, Si and Mg) are relevant. 
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Fig. 3. — Upper panel : residual continuum CLV's 
at 3499.47 A derived from the 1-D (dashed) and 
the 3-D model (solid) with varied Fe abundances 
of ±0.3 dex. The unaltered data are highlighted. 
Lower panel : Difference of the data in the upper 
panel from the calculation with epo = ±0.3 dex. 

In order to investigate the impact of changes 
of the opacity on the Center-to-Limb variation we 
have calculated the R-CLV's for the 3-D and our 
reference 1-D models at 3499. 47 A with two differ- 
ent Fe abundances (±0.3 dex). The purpose of the 
test is to demonstrate the general effect of opacity 
variations that can come from different sources, 



i.e. uncertainties of abundances and uncertainties 
of bound-free cross-sections of all relevant species 
(not only iron). However, to simplify the proce- 
dure we have modified only the abundance of iron 
which stands for the cumulative effect of all un- 
certainties. In the example the total opacity is 
increased by 50% and decreased by 22%, respec- 
tively. 

Increased opacity, i.e. increased iron abun- 
dances, results in large negative residuals while 
decreased opacity produces large positive residu- 
als (Fig. [3|). Both models are affected in a simi- 
lar way, but the strength of the effect is slightly 
smaller for the 3-D calculation by about 20% (cf. 
Fig. [3J lower panel). A change in opacity has a 
significant effect on the CLV but it does not elim- 
inate the discrepancies. 

To estimate the effect of a varied temperature 
gradient on the R-CLV's we have calculated the 
CLV at 3499.47 A for two artificially modified 1- 
D models (Fig. [4]). The temperature structure 
around tr oss = 2/3 is changed such that the gra- 
dient in temperature is increased and decreased 
by 1%, respectively. At fj. = 0.2, i.e. the po- 
sition of the largest discrepancy, the residual of 
0.0085 is changed by « 0.0035, i.e. by roughly 
1/3, indicating a maximum error of the 1-D tem- 
perature gradient of about 3%. Again, a simple 
change does not lead to perfect agreement, espe- 
cially when more than one frequency is considered. 
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Fig. 4.— Residual continuum CLV's at 3499. 47 A 
from the 1-D model with modified temperature 
fields. The temperature gradient at t « 2/3 is 
changed by +1% (upper curve) and -1% (lower 
curve). The unaltered data are highlighted. 

Finally, we compare the CLV in the continuum 
for the average ('horizontal' and over time) 3-D 
model with the exact data derived from the ra- 
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Fig. 5. — Residual CLV's in the continuum, i.e. the difference of the observed and the synthetic normalized 
CLV's for the 3-D model (solid) and its 1-D approximation (dashed). In comparison, the approximation 
in 1-D, i.e. the 'horizontal' time average of the 3-D model, fails to reproduce the observed Center-to-Limb 
variation at all. The maximum absolute residual is up to 5 times larger, i.e. -0.15 (1-D) vs. 0.03 (3-D). 



diation transfer in 3-D. Fig. shows the residual 
CLV's for both models. The discrepancies with 
the observations are much more severe for the av- 
erage 3-D model and it becomes obvious that it 
does not represent the original 3-D time series at 
all. Although a 1-D representation would obvi- 
ously be highly desirable because it would allow 
to quickly calculate spectra by means of a 1-D ra- 
diation transfer code, this turns out to be a very 
poor approximation in this case. 



Avres et al.l ((2006) have carefully investigated 



the rotational- vibrational bands of carbon monox- 
ide (CO) in the solar spectrum and have de- 
rived oxygen abun dances from three mo dels, i.e. 
the Fal C model (|Fontenla et alJ Il993l ). a 1-D 



model that is especially adapted to match the 
Center-to-Limb variation of the CO-bands (CO- 
mosphere), and from the averaged 3-D time se- 
ries. In all three cases, temperature fluctuations 
are accounted for in a so-called 1.5-D approxi- 
mation, in which profiles from 5 different tem- 
perature structu res are averaged. B y assuming a 
C/O ratio of 0.5. lAvres et"aH (|2006l ) derive a high 
oxygen abundance close to the "old" value from 
iGrevesse fc Sauvall (|l998t ) from both the Fal C and 
the COmosphere model, discarding the low oxygen 
abundance derived from the mean 3-D model be- 
cause its temperature gradient is too steep around 
To. 5 Mm ~ 1 an d fails to reproduce the observed 



Center-to-Limb variations. 

Our current study documents that the mean 3- 
D model is not a valid approximation of the 3-D 
time series, and therefore its performance cannot 
be taken as indicative of the performance of the 
3-D model, and in particular of its temperature 
profile. We find that the Center-to-Limb variation 
of the continuum predicted by the 3-D simulation 
matches reasonably well (i.e. similar to the best 
1-D model in our study) the observations. The re- 
sults by Scott et al., based on 3-D radiative trans- 
fer on the same hydrodynamical simulations used 
here, indicate that the observed CO ro- vibrational 
lines are consistent with the low oxygen and car- 
bon abundances. Our results show that there is 
no reason to distrust the 3-D-based abundances 
on the basis of the simulations having a wrong 
thermal profile. 

3.2. Lines 

We study the Center-to-Limb variation of a 
number of lines by comparing observations of the 
quiet Sun taken at 6 different heliocentric angles to 
synthetic profiles derived from 3-D and 1-D mod- 
els. The observation s are described in detail by 
lAllende Prieto et al. ( 2004 ) and were previously 
used for the investigation of inelastic collisions 
with neutral hydrogen on oxygen lines formed un- 
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der non- LTE condition^. The observations cover 
8 spectral regions obtained at 6 different positions 
on the Sun. The first 5 slit positions are centered 
at heliocentric angles of \i = cos9 = 1.00, 0.97, 
0.87, 0.71 and 0.50. The last position varies be- 
tween [i = 0.26 and 0.17 for different wavelength 
regions. This translates to distances of the slit cen- 
ter from the limb of the Sun in arcmin of 16.00', 
12.11', 8.11', 4.74', 2.14', 0.54' and 0.24', assuming 
a diameter of the Sun of 31.99'. For both of these 
last positions the slit extends beyond the solar disk 
and the center of the illuminated slit corresponds 
to \i = 0.34 and 0.31 (0.96' and 0.78'). 

We have calculated a variety of line profiles 
for the 6 positions defined by the center (in /x) 
of the illuminated slit. Although the slit length, 
160arcsec, is rather large, test calculations show 
that averaging the spectrum from six discrete \x- 
angles spanning the slit length gives virtually the 
same equivalent width than the spectrum from the 
central /x. For /i = 0.5, the second last angle, the 
difference amounts to a marginal change of the 
\og-gf value of about 0.01. To further reduce the 
computational burden we have derived the aver- 
age 3-D profiles from calculations taking only 50 
(every other) of the 99 snapshots into account. 

We have selected 10 seemingly unblended lines 
from 5 different neutral ions. The list of lines is 
compiled in Table [TJ The log-gf values for most 
lines were adopt ed from laboratory me asurements 
at Oxford (e.g. iBlackwell et all (Il995f) and refer- 
ences therein) and by lO 'Brian et al. nfl99ll) . 

We are interested in how synthetic lines profiles 
deviate from observations as a function of the posi- 
tion angle \i for two reasons. First of all, any clear 
trend with \i would reveal shortcomings of the the- 
oretical model atmospheres similar to our findings 
presented in Sect. 13.11 But, arguably, even more 
relevant is the fact that any significant deviation 
(scatter) would add to the error bar attached to a 
line-based abundance determination. 

In our present study we compare synthetic line 
profiles from 3-D and 1-D models with the obser- 
vations. Due to the inherent deficiencies of the lat- 
ter models, i.e. no velocity fields and correspond- 
ingly narrow and symmetric line profiles, etc., we 
focus on line strengths and compare observed and 
synthetic line equivalent widths, rather than com- 
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Fig. 6. — Iron lines under consideration. Ob- 
served (dashed) and synthetic (solid) profiles are 
shown for the disk center. The grey areas mark 
the velocity ranges used for the calculation of the 
line equivalent widths. The zero of the velocity 
scale refers always to the center of the observed 
line profile, as approximately determined by poly- 
nomial fitting. The log-<?/ values are modified to 
match the observed line equivalent widths. For 
most lines the profiles match well. However, the 
synthetic profile of the line at 5247.0 A seems to 
be broader than the observed profile. The line at 
6200.3 A is marginally blended around +5kms _1 . 
The line at 6170.5 A is noticeable blended around 
+5kms _1 . For both lines the wavelength interval 
is decreased accordingly. 



5 Data available at http://hebe.as.utexas.edu/izana 
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Tabic 1: Lines 



Ion A 


It 


max(0) 


log-fl 1 / 


lOg 1 R, a d 


log 1 Stark 


log 1 vdW 


r 1 1 

[A] 


tn\ 
\" ) 


[deg] 




I'M 


(r) 


{(!■) 


Fel 5242.5 


56000 


75 


-0.970 


7.76 


-6.33 


-7.58 


Fel 5243.8 


56000 


75 


-1.050 


8.32 


-4.61 


-7.22 


Fel 5247.0 


56000 


75 


-4.946 


3.89 


-6.33 


-7.82 


Fel( e ) 6170.5 


77000 


80 


-0.380 


8.24 


-5.59 


-7.12 


FelW 6200.3 


206000 


75 


-2.437 


8.01 


-6.11 


-7.59 


Fel 7583.8 


176000 


80 


-1.880 


8.01 


-6.33 


-7.57 


CrI 5247.6 


56000 


75 


-1.627 


7.72 


-6.12 


-7.62 


Ni I 5248.4 


56000 


75 


-2.426 


7.92 


-4.64 


-7.76 


Si I 6125.0 


77000 


80 


-0.930 


... (9) 






Til 6126.2 


77000 


80 


-1.425 


6.85 


-6.35 


-7.73 



a cf. lAllende Prieto et al. | l|2004l). R' is the reso lving power measured relative to the FTS (-Rfts ~ 400 000) spectrum at the 
center disk provided bv lBrault fc Neckell J1987T) . 

"r = 7, where 7 is the damping constant (FWHM of a Lorentzian profile, see, e.g. Eq. (11.13) in Gray 1992), in rad s _1 . 
c r = 7/Ne, where Ne indicates the number density of electrons at a temperature of 10,000 K (cgs units). 
d T = 7/Nfj, where Nh is the hydrogen number density at a temperature of 10,000 K. 
e noticeably blend at +5kms _1 
* marginal blend at +5kms — 1 

9 Approximate values were adopted for this line, see, e.g. Gray 1992 



paring the line profiles in detail. To be able to 
detect weak deviations, we have devised the fol- 
lowing strategy. We have identified wavelength 
intervals around each line under consideration for 
the contribution to the line equivalent widths and 
have calculated series of synthetic line profiles in 
1-D and 3-D with varied log-<?/ values that encir- 
cle the observations with respect to their equiv- 
alent widths. That allowed us to determine by 
interpolation the log-gf value required to match 
the observed line equivalent widths separately for 
each position angle ("Best-Fit"). To keep interpo- 
lation errors at a marginal level we have applied 
a small step of A(log-g/) = 0.05 for these series 
of calculations. A simple normalization scheme 
has been applied. All profiles have been divided 
by the maximum intensity found in the vicinity 
of the line center(within ±15kms _1 ). We con- 
volved the synthetic profiles with a Gaussian as to 
mimic the instrumental profile (see Table 1). An 
additional Gaussian broadening is applied to the 
line profiles from the 1-D calculation to account 
for macro-turbulence; this value was adjusted for 
each line in order to reproduce the line profiles 
observed at the disk center. 



Finally we have translated variations of line 
strength into variations of abundance, i.e. we have 
identified A log-(g/) = Alog-e. This approxima- 
tion is valid because the impact of slight changes 
in a metal abundance on the continuum in the op- 
tical is marginal. Note that it is not the intent 
of this study to derive metal abundances from in- 
dividual lines. Such an endeavor would require a 
more careful consideration regarding line blends, 
continuum normalization, and non-LTE effects. 

All calculations described in this section are 
single-line calculations, i.e. no blends with atomic 
or molecular lines are accounted for. The obser- 
vations did not have information on the absolute 
wavelength scale (see Allende Prieto et al. 2004), 
but that is not important for our purposes and 
the velocity scales in Figs. [6] and [7] are rela- 
tive to center of the line profiles. The individ- 
ual synthetic profiles were convolved with a Gaus- 
sian profile to match the observed profiles (cf. Ta- 
ble [I]). We were generally able to achieve a better 
fit of the observations when slightly less broaden- 
ing was applied to the 3-D profiles (0.3% in case 
of Fel 5242.5). Since we know from previous in- 
vestigations that the theoretical profiles derived 
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Fig. 7. — Non-iron lines under consideration. Ob- 
served (dashed) and synthetic (solid) profiles are 
shown for the center disk. The grey areas mark the 
wavelength ranges used for the calculation of the 
line equivalent widths. The \og-gf values are mod- 
ified to match the observed line equivalent widths. 
For most lines the profiles match well. However, 
the synthetic profile of the Sil line at 6125.0 A is 
noticeably narrower than observed. 



from 3-D Hydro-models match the observations 
well, we argue that the resolution of the obser- 
vat ions is actually slig htly higher than estimated 
by lAllende Prieto et al. I (|2004h . An alternative 
explanation would be that the amplitude of the 
velocity field in the models is too high. Such a 
finding, if confirmed, deserves a deeper investiga- 
tion but is beyond the scope of this study since 
line equivalent widths are only marginally (if at 
all) affected. 

We introduce the lines under consideration by 
showing the observed center-disk line profiles and 
the "Best-Fits" derived from the 3-D calculations 
of the six Fel lines and the four lines from other 
ions in Figs. [5] and respectively. In Fig. [5] we 
exemplify the fitting process by means of the Fe I 
line at 5242.5 A and show the relative difference 
between the observation and a variety of model 
calculations for all 6 angles under consideration. 
The "Best-Fit" \og-gf values are derived by inter- 
polation to match the observed equivalent widths 
from the spectral region around the line profile. 

We have obtained "Best-Fits" for all 10 lines 
(cf. Table [1} and present the log-gf values as a 
function of /i in Fig. [51 Be reminded that the 
aim of this study is not the measurement of ab- 
solute abundances: we focus on relative numbers 
and normalize our results with respect to the disk 
center (^ = 1). 

For improved readability we subdivide our find- 
ings presented in Fig. [9] into 4 distinct groups, i.e. 
iron/non-iron lines and 1-D/3-D calculations, re- 
spectively. We focus our discussion on the first 
five data points because we have some indications 
that the data obtained for the shallowest angle is 
less trustworthy than the data from the other an- 
gles: i) the relative contribution of scattered light 
was estimated from the comparison of the center- 
disk spectrum w i th the FTS spectrum taken from 
Brault fe Neckell ( 1987 ). and the outer- most posi- 



tion was the only one for which the entire slit was 
not illuminated, ii) for all 10 lines the fit of the line 
profiles for this particular angle is the worst (cf. 
Fig. [8]) and Hi) the scatter in our data presented 
in Fig. [9] is the largest for this angle. Fortunately, 
the flux integration is naturally biased towards the 
center of the disk. 

We find this systematic behavior for all 6 iron 
lines: Alog(e) is larger or equal in 1-D compared 
to 3-D, for all but one line (Fel 6170.5 A) Alog(e) 
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Fig. 8.— Fit of the Fel line at 5242.5 A with syn- 
thetic profiles derived from the 3-D Hydro model. 
Upper panel: Normalized profiles of the center 
disk observation (dashed) and the "Best-Fit" syn- 
thetic profile (solid). Lower panels: Relative dif- 
ference of the observation and three synthetic pro- 
files (dashed), i.e. three different log ^/-values for 
the 6 angles under consideration. The equivalent 
widths of the "Best-Fits" (solid) match the ob- 
served equivalent widths. The wavelength range 
considered for the equivalent width is highlighted 
in grey. The fit improves when a slightly higher 
resolution (by « 0.3%) is assumed. However, line 
equivalent widths are only marginally affected. 



is positive or zero for the 1-D calculations, and 
Alog(e) is negative or zero for all 3-D calculations. 
The Fel line at 6170.5 A stands out in both com- 
parisons. In 1-D it is the only line with a nega- 
tive Alog(e) and in 3-D it shows the by far largest 
negative Alog(e). This might be related to the 
noticeable line blend (cf. Fig.[6|). 

The iron lines calculated in 3-D indicate a uni- 
form trend of decreased \og-gf values with in- 
creased distance from the center-disk. The average 
decrease at /i = 0.5 for this group is -0.015 (Fel 
6170.5 A excluded). From the 1-D calculations we 
derive the opposite trend for the same group of 
lines and obtain an average of 0.103. Obviously, 
the 3-D model performs significantly better than 
the 1-D reference model regarding the center-to- 
limb variation of Fe I lines, even when equivalent 
widths, and not line asymmetries or shifts, are 
considered. 

For these five Fel lines we obtain an average 
difference (1-D vs. 3-D) of 0.12 at fj, = 0.5. To 
estimate the impact on abundance determinations 
based on solar fluxes we apply a 3-point Gaus- 
sian integration, neglecting the shallowest angle 
at fi = 0.11 (which has, by far, the smallest in- 
tegration weight) for which we have no data, and 
assuming that the good agreement between the 
1-D and the 3-D calculations for the central ray 
implies an equally good agreement for the first 
angle at \x = 0.89. These estimates lead to an 
abundance correction of approximately 0.06 dex 
between 1-D and 3-D mode ls due to their differ- 
ent center-to- limb variation. Asplund et alJ ( 200Ch 
found a similar correction from the comparison of 
1-D and 3-D line profiles at the disk center. 

For the 4 non-iron lines we find a uniform trend 
of increasing log-gf values with decreasing /i for 
both, the 1-D and the 3-D dataset. The systematic 
behavior is similar to what we find for the iron 
lines, but now the performance of the 1-D and 3- 
D models is similar, and the offsets are in the same 
sense: larger abundances would be found towards 
the limb for both models. 

4. Conclusion 

The photosphere of cool stars and the Sun 
can be described by stellar atmospheres in 1-D 
and 3-D. Since the 3-D models add more realis- 
tic physics, i.e. the hydrodynamic description of 



14 



0.3 
0.2 

0.1 

0.0 

-0.1 

-0.2 

0.3 
0.2 

0.1 

0.0 

-0.1 

^-0.2 

j? 0.3 
< 0.2 

0.1 

0.0 

-0.1 

-0.2 

0.3 
0.2 

0.1 

0.0 

-0.1 

-0.2 

-0.3 



0.0 



Fe I 5247.0 
Fe I 5242.5 
Fe I 7583.8 

Fe I 6200.3 
Fe I 5243.8 



1-D 



=5 _ _ 



Fe I 6170.5 



3-D 



Fe I 5243.8 
Fe I 5247.0 



Fe I 5242.5 
Fe I 6200.3 
Fe I 7583.8 

Fe I 6170.5 




1-D 



Cr I 5247.6 
Ni I 5248.4 
Ti I 6126.2 



Si I 6125.0 



Cr I 5247.6 ^ 
Ni I 5248.4 
~ Si I 6125.0 -== 




3-DE 


= Til 6126.2 






— A 


A 


A A AE 


= • □ 


O 


□ o E 


E 0.031 0.106 

-iii 


0.222 


0.394 0.247 : 

iii- 



0.5 



1.0 



Fig. 9. — Relative abundance variations with respect to the center disk (/z=l) from the 1-D (dashed) and 
the 3-D calculations, respectively. The upper two panels show the six Fel lines, the lower two panels show 
the four lines from the other elements. Line designations are given in the left part of the plot. In the second 
panel from the top, the curves for the iron lines at 5242.5 and 6200.3 A overlap with each other for the two 
positions closest to the limb. The bottom panel also indicates the six ^-angles of the observations (triangles, 
/i=l not shown, shallowest angles varies slightly) and the /x-angles used in a 2-point (squares) and a 3-point 
(circles) Gaussian flux integration and their respective integration weights (J2 w ^ = 0.5). 
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the gas, they can be seen truly as an advancement 
over the 1-D models. However, this refinement in- 
creases the computational effort by many orders 
of magnitude. In fact, the computational work- 
load becomes so demanding, that the description 
of the radiation field has to be cut back to very few 
frequencies, i.e. to a rudimentary level that had 
been surpassed by 1-D models over 30 years ago. 
Overall we are left with the astonishing situation 
that a stellar photosphere can be modeled by ei- 
ther an accurate description of the radiation field 
with the help of a makeshift account of stellar con- 
vection (Mixing- length theory), or by an accurate 
description of the hydrodynamic properties aug- 
mented by a rudimentary account of the radiation 
field. 

It is evident that individual line profiles can be 
described to a much higher degree and without 
any artificial micro- or macro-turbulcncc by the 
3-D Hydro models, as the simulations account for 
Doppler-shifts from differential motions within the 
atmosphere. We know from detailed investigations 
of line profiles that the velocity field is described 
quite accurately and that the residuals of the fit- 
tings to line profiles are reduced by about a factor 
of 10. However, it is not obvious how the 3-D 
models compare to their 1-D counterparts when it 
comes to reproduce spectral energy distributions 
and line strengths. 

We study the solar Center-to-Limb variation for 
several lines and continua, to probe the tempera- 
ture structure of 3-D models. The work is facili- 
tated by the new code Asset, which allows for the 
fast and accurate calculation of spectra from 3-D 
structures. In compar i sons to other progr a ms (e .g. 
Asplund et all (|2000h ; iLudwig fc Steffenl (120071 )). 
the attributes of the new code are a greater 
versatility, i.e. the ability to handle arbitrarily 
complicated lines blends on top of non-constant 
background opacities, higher accuracy due to the 
proper incorporation of scattering and improved 
(higher-order) interpolation schemes, and a higher 
computational speed. 

In our study we find that regarding center-to- 
limb variations, the overall shortcomings of the 3- 
D model are roughly comparable to the shortcom- 
ings of the 1-D models. Firstly, we conclude from 
the investigation of the continuum layers that the 
models' temperature gradient is too steep around 
r « 2/3. This behavior is more pronounced for 



the 3-D model which shows a drop in intensity 
(with (i) that is about twice the size of the drop 
displayed by our reference 1-D model, but at the 
same time smaller than the discrepancies found for 
two other (newer!) 1-D structures. Secondly, the 
line profiles for different position angles on the Sun 
cannot be reproduced by a single abundance. For 
Fe I lines, the abundance variation between the 
disk center and \x = 0.5 is about 0.1 dex for our 
reference 1-D model, but only 0.015 dex (and with 
the opposite sign) for the 3-D simulations, albeit 
the calculations for lines of other neutral species 
suggest a more balanced outcome. 

Overall we conclude that the 1-D and the 3-D 
models match the observed temperature structure 
to a similar degree of accuracy. This is somewhat 
surprising but it might be that the improved de- 
scription of the convective energy transport is off- 
set by deficiencies introduced by the poor radia- 
tion transfer. Once new Hydro models based on an 
upgraded radiation transfer scheme (i.e. more fre- 
quencies and angles, better frequency binning) be- 
come available in the near future (Asplund, priv. 
comm.), we will be able to test this hypothesis. It 
will become clear whether focusing on refining the 
radiation transfer will be enough to achieve better 
agreement with observations, or the hydrodynam- 
ics needs to be improved as well. 
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